Melting at dislocations and grain boundaries: A Phase Field Crystal study 
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Dislocation and grain boundary melting are studied in three dimensions using the Phase Field 
Crystal method. Isolated dislocations are found to melt radially outward from their core, as the 
localized excess elastic energy drives a power law divergence in the melt radius. Dislocations within 
low-to-mid angle grain boundaries melt similarly until an angle-dependent first order wetting tran- 
sition occurs when neighboring melted regions coalesce. High angle boundaries are treated within a 
screening approximation, and issues related to ensembles, metastability, and grain size are discussed. 
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Freezing and melting transitions do not exhibit the 
range of universal behavior associated with continuous 
phase transitions and largely for this reason have eluded 
a unified theoretical description. The nature of a given 
melting transition may depend sensitively on the details 
of the system and experiment, and can involve many dis- 
tinct processes both within and between multiple forms 
of excitations. For example, melting may occur abruptly 
and discontinuously at the melting temperature T^, or 
it may initiate well below Tm at surfaces and/or internal 
defects and proceed up to T^. This process of premelting 
has been studied extensively for surfaces P, |3l and is rel- 
atively well understood, but limited and inconsistent ex- 
perimental evidence for melting at dislocations and grain 
boundaries leaves a number of issues unresolved. 

A recent study of colloidal crystals has verified that 
premelting does occur at vacancies, dislocations, and 
grain boundaries, and has provided measurements of the 
localized premelting behavior below Tm [Sj. The con- 
ditions which determine whether premelting will occur 
continuously or discontinuously, and whether the width 
of the premelted region diverges are not fully understood. 
Grain boundaries in Al have been found to liquify very 
near Tm, and the width of the melted layer appears to 
diverge |^]. Discontinuous jumps in grain boundary dif- 
fusion coefficients mobility 0, and shear resistance 
Q have been found in other metals. 

Theoretical studies have been based on either ex- 
plicitly atomistic methods such as Molecular Dynamics 
[9|, llO| and Monte Carlo or on continuum phase 

field models with uniform phases [13, [H, In t^iis 

study, dislocation and grain boundary melting are ex- 
amined using the Phase Field Crystal (PFC) method 
[l5| , which extends the phase field approach to the level 
of atomistic resolution. This permits straightforward 
identification of stable equilibrium and metastable non- 
equilibrium atomic structures, while inherently including 
crystal symmetry and orientation, elasticity/plasticity, 
and the individual dislocations which compose the grain 
boundaries. Our description will be most applicable to 



hard-sphere/colloidal systems and possibly simple met- 
als. 

The melting behavior of dislocation pairs and symmet- 
ric tilt grain boundaries of 61 = 4°, 8°, 16°, 24°, 32°, and 
44° are examined numerically for a simple PFC model 
with bcc symmetry. Analytic results are derived for iso- 
lated dislocations and low boundaries by combining 
the PFC equations with continuum linear elasticity. A 
screening approximation is outlined for high angle bound- 
aries, though somewhat surprisingly, the low 9 descrip- 
tion is found to remain reasonably accurate for high 9. 

As shown in reference the PFC free energy can 
be derived from the Ramakrishnan-Yussouff free energy 
functional of classical density functional theory IT], [l^ . 
Here we give the final form: 
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where F = {T — To)/ {pkBTL'^), To is the free energy 
functional at constant density, p is the average num- 
ber density, ks is Boltzmann's constant, T is temper- 
ature, L = \j 21(741/(72, Ci is the i point direct correla- 
tion function of the reference hquid state, = 1 — pCo, 



= piC2y/4\C4\, 



x/L and v accounts for the 



lowest order contribution from C3. n — (p ~ Pl)/pl is 
the scaled number density, where p is the local density 
variable and p^ — p is the liquid number density. Clas- 
sical density functional theory has been used to examine 
surface melting [J [l^l , but not grain boundary melting, 
presumably due to the complexity of the solid-solid inter- 
face and the more demanding system size requirements. 
The dynamics are given in dimensionless form by 



dn 
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dr Sn 

where {r]{fi,Ti)ri{r2,T2)) = MV^S{fi ~r2)S{Ti -T2) and 
M = 3{i5L'^)^. This form imposes a constant density and 
is consistent with the canonical ensemble. 

A semi-implicit pseudospectral algorithm was used to 
solve Eq. ([2|) in three dimensions with periodic boundary 



FIG. 1: Numerically measured local melt radius Rm around 
an edge dislocation in a bcc crystal as a function of temper- 
ature (units of lattice constant a, AB < 0.025 values ob- 
tained by extrapolation of 4>if))- Ro ~ 0.2812a is the zero 
temperature radius, determined by best fit. Inset: cross- 
sectional images of n{x, y, z) from simulations at AB — 
0.02165, 0.02511, 0.02569, 0.02627, 0.02656, 0.02685 from a) to 
f), showing melting at a dislocation core. 



conditions. The parameters used were Ax — 0.976031, 
At = 0.5, ^ V3/3, v ^ 3^^^/2, and M = 0.002. A 
system size V — (512Aa;)'^ = (56a)'^ was used for the 
dislocation pair and 4° grain boundary pair, while V = 
(256 Ax) = (28a) was used for all other grain boundary 
pairs, where a = 8.9237 is the bcc lattice constant. Finite 
size effects increase as 9 decreases, but were found to 
be small for all grain boundaries studied. The gaussian 
width or mean square displacement (D) of each localized 
density peak was monitored as the temperature AB = 
B^ — i?" was increased toward the melting point. The 
local crystallinity (f> has been defined as (/)(r) = {D{f) — 
Dx)/{Dl — Dx), where Dx is the equilibrium D of the 
crystal phase and I?l is that of the liquid phase, (f) = 1/2 
specifies a liquid-solid interface. 

The radius of melted region around a dislocation core 
Rm was first measured in this manner for an edge dis- 
location pair as the temperature was raised toward the 
bulk melting temperature AB„i- The results are shown 
in Fig. [1] where the data is plotted as (i?,„ -I- i?o)~^ 
vs AB to demonstrate that is consistent with a 
{AB„i — AB)^^/^ form which will be derived later. i?o 
is an offset related to the finite size of the dislocation 
core at zero temperature. The fit to this form predicts 
a bulk melting temperature ABm — 0.0270 which is 
in good agreement with the directly measured value of 
ABm — 0.0278. The upper inset of Fig. [1] shows melting 
around an edge dislocation as AB ABm- 

Figure [2] shows the progression of melting at 8° and 
44° grain boundaries. Low angle boundaries were found 
to first melt radially at each dislocation core, until the 
melted regions of neighboring dislocations coalesce and a 
uniform wetting layer is formed along the boundary. Indi- 




FIG. 2: Laterally averaged crystallinity parameter (j) and 
cross-sections of n{x, y, z) for 8° (upper) and 44° (lower) grain 
boundaries at AB = 0.02165,0.02627,0.02656,0.02685 and 
AB = 0.02425,0.02468,0.02540,0.02656 from left to right, 
respectively. 

vidual dislocations cannot be distinguished in high angle 
boundaries, and melting in this case was found to pro- 
ceed by uniform disordering along the boundary rather 
than by local radial melting. Interfacial roughening due 
to thermal fluctuations is negligible in all simulations pre- 
sented here. 

The dependence of the width of the wetting layer (or 
the liquid volume fraction of the system) on AB is shown 
in the inset of Fig. [3] for various grain boundary angles. 
In all cases the width remains narrow and the bound- 
ary relatively dry until above the solidus, at which point 
a discontinuous jump is observed at some characteristic 
wetting temperature AB^et- The dependence of AB^et 
on 9 is shown in the main plot of Fig. [3l The fit lines 
will be discussed in the following, though the axes of 
the plot reveal already that our predicted form will be 
AB^et oc sin^ 9. 

Based on these simulation results, we have developed 
a theory of dislocation-driven melting, which is easily ex- 
tended to low angle grain boundaries. The low angle 
results are shown to remain accurate for all but the high- 
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FIG. 3: Grain boundary wetting temperature vs. sin Q. Fit 
lines are discussed in the text. Inset: Liquid volume fraction 
vs AB for a various grain boundaries. The dotted line cor- 
responds to the liquid volume fraction predicted for the 8° 
boundary based on Eq. ([SJ only (ignoring coexistence). 



est 9 where the dislocation spacing d reaches the order of 
the burger's vector h. A screening approximation for the 
spatial grain boundary energy is found to be more ap- 
plicable for very large 0, with a gradual crossover taking 
place between these two regimes. 

According to continuum elasticity [i^l , the radially av- 
eraged elastic energy density per length of dislocation line 
in a three-dimensional isotropic solid is F^i — a^/R^, 
where /x is the shear modulus. For a screw dislocation 
as = 6^/47r^ and for an edge dislocation = as/(l — ct) 
where tr is the Poisson ratio. If we assume this result 
to hold for a dislocation in the PFC model, at distances 
approaching the core region, then can be directly cal- 
culated by determining the distance at which F^i is suffi- 
ciently large to destabilize the crystalline phase, melting 
the dislocation core. 

It will be assumed that n can be represented in a 
one mode approximation for a bcc lattice, i.e., n{r) — 
A (cos qx cos qy + cos qx cos qz -\- cos qy cos qz). Substi- 
tuting n into Eq. ([!]) and minimizing with respect to 
q gives, 
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where Af-^ = (F ~~ F]^)/V, Fj^ is the free energy of the 
liquid, V = (27r/g)^, and qmin = •\/2/2. The shear mod- 
ulus can be estimated in the one mode approximation 
by setting n(a;, y, z) — > n{x + Qy, y, z) and expanding F 
in C such that F = F(C ^ + ^JiC,'^ ^ . This proce- 
dure gives ^/ksTL'^p — A'^B'^/S. The total (dimension- 
less) free energy of the system with a dislocation is then 
A/^ -I- Fei/ksTpL'', which can be written 

Af=-(AB + ^]A^--A' + —A\ (4) 
8 I i?2 / 8 256 ^ ^ 



where R = R/b and Es = B'/{12tt'^), = Es/{1 - a) 
for screw and edge dislocations respectively. Equation 
dH) indicates that the elastic energy 'shifts' the effective 
temperature AB by an amount E/R?. The implication 
is that the liquid-solid transition is shifted and instead 
of occurring when AB = Ai?,„ occurs when when AB -|- 
E/R^ = AB„i. Thus the premelt radius can be written 



i?™ = ^E/{AB„, - AB) 



(5) 



or = (AB„, - AB)IE. 

As shown in Fig. [H this form is consistent with the 
simulation results for edge dislocation pairs, though the 
predicted slope {—l/Eg) is in error by a factor of nearly 
five. A more definitive test would require additional data 
very near ABm, a region increasingly difficult to access 
due to system size requirements. 

Extending this result to low angle boundaries, we can 
estimate the grain boundary wetting temperature AB^et 
where neighboring dislocations coalesce by setting Rm — 
d/2 = 1/(2 sing). Substituting Eq. ^ for R^ gives 



AB„ 



ABr, 



AEsin-'e 



(6) 



which is in good agreement with the data shown in Fig. [S] 
This result should lose accuracy as 9 increases since the 
dislocation energies gradually deviate from the isolated 
dislocation result as 6* ^ 9max- The agreement up to 
9 ~ 32° is somewhat unexpected as the superposition 
generally loses accuracy for 9 > 10°. The best fit line 
predicts AB^ = 0.0272, again near the measured value. 

The solid line in Fig. [3] corresponds to the fit line from 
Fig. m set equal to 1/(2 sin9) and solved for AB^et- The 
agreement here clearly indicates that the wetting of low 
and mid-angle grain boundaries is accurately described 
by the coalesence of radially melted regions around nearly 
isolated dislocations. 

In the limit of large 9 {d — > 0), the grain bound- 
ary energy becomes increasingly uniform along its length 
(see Fig. ^ and can no longer be described in terms 
of individual dislocations. We expect that elastic fields 
at long distances are screened by closely spaced dis- 
locations, giving rise to exponentially decaying spatial 
grain boundary energy. Indeed, direct analysis of free 
energy data from simulations indicates that such an 
exponential form is qualitatively correct. Solving for 
ABu,et using Fei oc e"^/'' rather than F^i cx 1/i?^ gives 
ABm — AB^et oc e"-'''^^^ This is the form of the wide 
dashed line in Fig. [3l which appears to be more appro- 
priate for large 9. 

Some comments concerning the influence of liquid-solid 
coexistence and the canonical ensemble (i.e., conserved 
density) on grain boundary melting may be helpful at 
this point. The equilibrium state for a simple system 
with a grain boundary is most generally either dry if 
Fgb < 2Fu+£{Fl-Fx) or wet ifFgt > 2i^j,+^(Fi-Fx), 
where Fgb is the grain boundary energy, Fis is the energy 
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of a liquid-solid interface, and I is the width of the liq- 
uid region in the wet state. If the wet state becomes 
favorable below the melting temperature, then a grain 
boundary wetting transition occurs. In the canonical en- 
semble as examined here, the effects of liquid-solid coex- 
istence and the subsequent shifts in density of the two 
phases above the solidus ABsoi modify this heuristic ar- 
gument. Now ABm, the temperature at which Fl = Fx, 
is straddled by a coexistence region. As AB A_B,„ the 
system first encounters a solidus above which some vol- 
ume fraction of liquid will minimize the overall F. For 
the grain boundary pair geometry, the equilibrium state 
above ABsoi is one with a uniform volume of liquid occu- 
pying each boundary region. Therefore, an equilibrium 
first order wetting transition will occur at ABgoi as long 
as the grain size is not excessively small. Above ABgoi, 
the liquid layer width will be controlled by coexistence 
rather than local defect energies, since the elastic fields 
of the grains largely decouple {Fei 0) upon wetting. 

The results presented here show no wetting or strong 
premelting for AB < ABgoi, and the equilibrium wet- 
ting transition is not observed. Instead, a ^-dependent 
discontinuous transition from the metastable dry bound- 
ary state to the equilibrium wet state occurs at AB^oi < 
ABi^et < ABjn, as shown in the inset of Fig. [3l This is 
because the wetted state is not nucleated in observable 
times until i?„i has grown sufficiently large to coalesce 
and the free energy barrier approaches zero. The disloca- 
tions and/or grain boundaries act as nucleation sites for 
the liquid above ABsoi, creating well-defined nonequi- 
librium paths from the metastable dry state to the F 
minimizing wet state (which all must conserve p). 

The condition for wetting in the canonical ensemble 
involves the grain size Lg, such that wetting can be 
supressed to temperatures above AB^oi when Lg is fi- 
nite. The condition can be written Fgf, + LgAFx > 
2Fis+eAFc where AFx = Fx[p] - Fx[px] and AFc = 
FlIpl] - Fx[px]- Here £ = {px - p)/{px - Pl), P is the 
conserved system density, and px and pL are the shifted 
coexistence densities of the solid and liquid phases, re- 
spectively. For AB < ABsoi , if we assume that AFx = 
and PL — Px ~ P (this is not the case when premelting is 
strong below ABgoi), we recover the original inequality 
Fgb > 2Fis+£{Fl — Fx) and Lg is not a significant factor. 
In the limit — s- oo, the wetting condition will always 
be satisfied for AB > ABsoi and the equilibrium transi- 
tion occurs at ABsoi- As Lg decreases, the equilibrium 
wetting transition is shifted to higher AB. 

Two-dimensional systems show qualitatively similar 
results to those presented above. Preliminary simula- 
tions have also been conducted in the grand canonical 
ensemble (i.e., non-conserved density), where the above 
complications with liquid-solid coexistence due to den- 



sity conservation can be avoided. Initial results indicate 
that there is no strong wetting transtion in this case, as 
the atoms near grain boundaries continuously delocalize 
toward ABm but do not liquify and decouple the grains 
until above ABm- 
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